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Abstract 

We investigate the geometry, stability, electronic structure and optical proper- 
ties of C24H12 coronenes encapsulated in a single-wall (19,0) carbon nanotube. By 
an adequate combination of advanced electronic-structure techniques, involving 
weak and van der Waals interaction, as well as many -body effects for establishing 
electronic properties and excitations, we have accurately characterized this hybrid 
carbon nanostructure, which arises as a promising candidate for opto-electronic 
nanodevices. In particular, we show that the structure of the stacked coronenes in- 
side the nanotube is characterized by a rotation of every coronene with respect to 
its neighbors through van der Waals interaction, which is of paramount importance 
in these systems. We also suggest a tentative modification of the system in order 
this particular rotation to be observed experimentally. A comparison between the 
calculated many-body excitation spectrum of the systems involved reveals a pro- 
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nounced optical red-shift with respect to the coronene-stacking gas-phase. The 
origin of this red-shift is explained in terms of the confinement of the coronene 
molecules inside the nanotube, showing an excellent agreement with the available 
experimental evidence. 

1. Introduction 

In recent years, substantial effort has been devoted to the development of new 
nanostructures involving standard molecules, which states as a key concept for 
their applicability in state-of-the-art future nanodevices [2-4]. This burgeoning 
researching field gives rise to the building of more and more complex struc- 
tures |51 [6l, which are often inexpensive, versatile and stable materials that can 
traditionally be used at low temperatures, and are nowadays becoming attractive 
for high-temperature regimes. On this basis, the self- assembling of molecules 
through confinement within carbon nanotubes E [H offers a wide variety of in- 
teresting supramolecular structures. These systems are constituted of building 
blocks whose stability is ensured by weak intermolecular interactions. Among 
these mentioned interactions, the most fundamental and dominant one is the van 
der Waals attraction (vdW) [0 [TOj . It allows the assembling of molecules - or 
groups of molecules - together without substantially affecting the basic proper- 
ties of each molecule individually. In fact, supramolecular chemistry 0TJ, as well 
as many biological systems [TT2l . are based on these assembling of building blocks 
through vdW interaction. Nevertheless, despite the fundamental relevance of such 
interaction, its treatment is still very complicated to handle, especially from a first- 
principles framework. Some methods have been proposed in literature on the basis 
of Density Functional Theory (DFT) [13-16], with relatively succesful results, but 
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which are often limited by the size of the considered systems, and the computa- 
tional resources. 

In particular, the case of self- assembled stacked coronene molecules inside 
carbon nanotubes is of particular interest since: i) first, the coronene molecules 
can be considered as small pieces of graphene, and as such, present interesting 
electronic properties; ii) second, as planar molecules, coronene molecules can be 
easily arranged to form huge ID-like supramolecular assembling. This property 
may be exploited for electronic transport purposes along the carbon nanotube axis; 
and iii) finally, as we shall revise in this work, the confinement of coronenes inside 
carbon nanotubes exhibits a non neglectable effect on their electronic structure. 
Consequently, these changes in the electronic structure may be utilized (and even 
tuned) by varying the size of the nanotube, in order to create new and efficient 
opto -electronic devices. 

In this communication, motivated by recent experimental results by Okazaki 
and coworkers (T), we consider the case of 7r-stacked coronenes in a Single Walled 
Carbon Nanotube (SWCNT). With this study, we aim at giving new and deep in- 
sights in the structure and energetics, as well as the electronic properties, of such 
system. This complex has been partially characterized theoretically through stan- 
dard DFT calculations in previous literature, but without considering the van der 
Waals interaction, which is of paramount importance for the intrinsic cohesion of 
this system. We shall show that a poor (or null) consideration of this challenging 
interaction leads to an erroneous description and characterization of this structure. 
On the other hand, from the electronic structure point of view, "many-body" tech- 
niques are dramatically necessary to properly describe the molecular gap, as well 
as the SWCNT band structure, - and, therefore, the excitation spectrum -. To our 
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knowledge, such an accurate framework has not been considered before. 

In the following, in section 2 we will present our calculation methodology. We 
will recall our intermolecular perturbation theory plus DFT to determine the struc- 
tural configuration of the coronenes inside the nanotube. Then, we will present 
a "many-body" corrected DFT approach to determine accurate band- structure di- 
agrams, as well as the excitation spectrum. In section 3 we will deeply discuss 
the structure obtained; and in section 4, we will provide the electronic structure 
results. All the previous results combined will be used to justify a very recent 
experimental evidence, by which an optical red-shift effect of the self-assembled 
C 24 H 12 coronenes stacked inside the (19,0) SWCNT with respect to its diluted 
clean-phase is observed [1J. On the basis of our findings, future experimental 
work is proposed. 

2. Method 

In order to shed some light on the interpretation of the mentioned optical red- 
shift evidence 01, the geometrical and electronic structure, as well as the optical 
photoexcitation properties, of all these systems have been deeply analyzed from 
an accurate theoretical point of view by the combination of advanced theoretical 
techniques: the van der Waals interaction for the geometry and energetics, and the 
"many-body" effects for the effective correction of the electronic structure. On 
this basis, we have performed a full state-of-the-art theoretical characterization of 
all the relevant systems involved in the present study. 

2.1. The LCAO-S 2 +vdW formalism 

As starting point, we first study the isolated coronene dimer, in order to prop- 
erly describe, subsequently, the confinement of stacked coronenes inside the nan- 

4 



TILTING 




Stacking distance 



Minimum energy structure 




3 3,5 4 4,5 

Stacking distance (A) 



Stacking distance 



ROTATION 




Minimum energy structure 




3,5 4 4,5 5 
Stacking distance (A) 



MINIMUM ENERGY 
STRUCTURE 



30° 





3,5 4 4,5 

Stacking distance (A) 



Figure 1: (Color online) (Top panel right) Cohesive energy of an isolated C24H12 coronene dimer 
as a function of the intermolecular stacking distance for several tilt angles calculated in the LCAO- 
5 2 +vdW formalism. (Middle panel right) Ibid for several roration angles with respect the axis 
dimer - note that rotation angles of 40, 50 and 60 degrees are respectively equivalent to rotation 
angles of 20, 10 and degrees -. (Bottom panel right) Ibid for several tilt angles, considering a 
rotation angle of or 30 degrees with respect to the dimer axis. Sketchs about tilt and rotation 
angles, and minimum energy structures are shown in left subpanels. 
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otube. All these systems are dominated by weak and van der Waals interactions, 
which will be accounted along the whole characterization process. For this pur- 
pose we have used the LCAO-S* 2 + vdW ifPTl [LSI . This approach is based on 
the Density Functional Theory (DFT) in combination with an intermolecular per- 
turbation theory to describe weak and van der Waals (vdW) interactions. Each 
interacting subsystem, namely in this case the coronene and the carbon nanotube, 
are calculated separately by DFT. The DFT computational scheme, as well as the 
theoretical foundations underlying our calculations - a very efficient DFT local- 
ized orbital molecular dynamics technique (Fireball) -, have been described in 
full detail elsewhere [19-22], and here we will only summarize the main points. 

We first analyze all the systems involved in this study by using a self-consistent 
version of the Harris-Foulkes LDA functional 11231 124| - instead of the traditional 
Kohn-Sham (KS) functional based on the electronic density -, where the KS po- 
tential is calculated by approximating the total charge by a superposition of spher- 
ical charges around each atom. The Fireball simulation package uses a local- 
ized optimized minimal basis set ll2~5ll . and the self-consistency is achieved over 
the occupation numbers through the Harris functional [26J. Besides, the LDA 
exchange-correlation energy is calculated using the efficient multicenter weighted 
exchange-correlation density approximation (McWEDA) [|20ll2T| . To these DFT 
calculations we add "weak interactions", which can be seen as two opposite con- 
tributions. The first one, named "weak chemical" interaction is due to the small 
overlaps between electronic densities of the interacting subsystems. Therefore, 
this energy can be determined as an expansion of the wavefunctions and operators 
with respect to these overlaps. This expansion is based on a development in S 2 
(since in the weak- interacting case the overlaps are really small) of the S -1 / 2 term 
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appearing in the Lowdin orthogonalization, which induces a shift of the occupied 
eigenenergies of each independent subsystem, and leads to a repulsion energy 
between them. This shift can be defined as follows: 

& E-n ^ ^ 2 l^nmTmn T nm S mn \ + — ^ ^ |*Snm| i^n £-m)) 

m rn 



n 



where S^m and T mn are the overlapping and hopping integrals between eigenvec- 
tors m and n in the two interacting subsytems, respectively; i.e. each coronene 
molecule, or the SWCNT and the molecule. The effect of the hopping matrix 
elements, T mn , can be calculated in a standard intermolecular second-order per- 
turbation theory as: 
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6 T e m = J2^^- (2) 

Thus, we obtain the following "weak chemical" contribution to the interaction 
energy: 

E weak chemical = 

2 (5 s e n + 5 T e n ) + 2 y~] (5 s e m + 5 T e m ), (3) 

n=occ. m=occ. 

where a multiplicative factor of 2 has been included to take into account the spin 
degeneracy and only filled states are considered. Typically, in graphitic or hydro- 
carbonated materials this contribution is in general repulsive. 

The second contribution is the pure van der Waals interaction, which finds its 
origin in charge fluctuations in each subsystems, arising from oscillating dipoles 
whose interaction gives the attractive part of the cohesive energy. This interac- 
tion is treated in the dipolar approximation, and added in perturbations to the total 



energy of the system. The balance of the two contributions gives the equilibrium 
configuration of the system. This formalism ifTTl [T8l has already provided excel- 
lent results in the study of a wide range of graphitic materials [27], or encapsulated 
molecules in nanotubes ||28ll . 

2.2. Band- structure diagrams 

Based on the previously obtained optimized vdW geometries, density func- 
tional theory (DFT) and "many-body" calculations - via the quasi-particle Green 
function screened interaction approximation (GW) ll29ll . combined in the next sec- 
tion to the Bethe-Salpeter equation (BSE) [30-33] for the calculation of the excita- 
tion spectrum - have been carried out on all the systems, in order to provide deep 
insights on their electronic structure (band-structure diagrams) and optical pho- 
toabsorption behavior (excitation spectrum). By accounting these "many-body" 
corrections, this study goes much deeper in the theoretical characterization of 
these systems than in any previous literature Q. 

DFT implementations used in practice have some well-known limitations due 
to the approximate description of exchange and correlation effects between the 
electrons, such as the known underestimation of the magnitude of the energy gap 
in semiconductors [|34l . A similar error could be expected for the gap between oc- 
cupied and unoccupied states (HOMO-LUMO) gap in finite clusters and nanopar- 
ticles [|35l . Possible errors in the gap and in the asymptotic behavior of the XC 
potential in a semiconductor bulk could dramatically affect the energies of the 
excited electronic states and, thus, the band structure, photoabsorption spectrum 
or electron-phonon coupling. For this reason, the DFT results are effectively cor- 
rected here with accurate electronic structure results obtained from "many-body" 
perturbation theory. To explore the influence of quasi-particle and electron-hole 
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interaction effects in the band-structure diagrams of the different systems, we have 
used the GW (Green function - screened interaction approach) 11291 . For this pur- 
pose: i) as a starting point, we took the previously obtained vdW geometries and 
lattice parameters for establishing the ground-state electronic structure of all sys- 
tems with the plane-wave PwScf code [38J, using for this step a Local Density 
(LDA) XC-functional ll39l - it is well-known that most-improving quasi-particle 
corrections are obtained from local density electronic ground-states (see 051 and 
references therein) -. Additionally, semicore-corrected norm-conserving pseu- 
dopotentials ll36l l37l were used to model the ion-electron interaction. The Kohn- 
Sham one-electron states were expanded in a series of plane-waves with an en- 
ergy cutoff of 500 eV, and the Brillouin zone of all systems was /c-sampled with 
(12 x 1 x 1) Monkhorst-Pack grids, guaranteeing a full convergence in energy and 
density. The self-consistent density was determined by iterative diagonalization 
of the Kohn-Sham Hamiltonian at k B T = 0.1 eV, using Pulay mixing of densities, 
and all total energies were extrapolated to k B T = eV. Spin-polarized calculation 
tests were also carried out to check the influence of the spin in these systems. 
Ground-state calculations for the vdW-based geometries converged fast, and re- 
vealed small forces on every atom ranging between 0.1 and 0.2 eV/A, which man- 
ifests an excellent transferability of geometries between simulation codes. This 
justifies their use in the more accurate electronic structure calculations. Addition- 
ally, spin-polarized calculations did not show either any significant variation with 
respect to the spin-unpolarized cases. Step ii) was performed with the Yambo 
code POl . and more than 200 unoccupied bands for each system were necessary 
to obtain self-consistent GW converged results. 
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2.3. Excitation spectrum 

A method extensively used to calculate the photoabsorption spectrum of nanos- 
tructures, surfaces and solids is the time-dependent density functional theory [|4T1 
l42l (TDDFT), which states a good performance for low electronically correlated 
systems through the time-evolution of the Kohn-Sham states. Nevertheless, as 
previously indicated, the implementations used in practice have some well-known 
limitations due to the approximate description of exchange and correlation effects 
between the electrons. For this reason, here we go beyond standard TDDFT by 
applying many -body perturbation theory via the quasi-particle approximation GW 
correction ll29l to the standard DFT, combined to the Bethe-Salpeter equation [30- 
33] for the calculation of the excitation spectrum, as follows: i) as a starting point, 
we consider the LGA+GW electronic quasi-particle corrections [29] to the elec- 
tronic structure for the ground-state of all the systems as calculated in the steps 

i) and ii) explained in the previous subsection; and ii) finally, the BSE [30-33] 
was solved for coupled electron-hole excitations ll33l . thereby accounting for the 
screened electron-hole attraction and the unscreened electron-hole exchange. Step 

ii) was also performed with the YAMBO code |40]|. This formalism, GW+BSE, 
includes, by construction, excitonic effects, which are more important in highly 
correlated systems, combining the excellent performance of the GW+LDA on the 
ground-state electronic structure calculations, with the high accuracy of the BSE 
for the calculation of excitations. 

3. vdW Geometries and Energetics 

The results of our study can be divided in two main aspects: i) first, we will 
consider the structural aspects of the C24H12 coronenes encapsulated in carbon 
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nanotubes due to the weak and van der Waals interactions. Coronene molecules 
stack into coaxial columns inside the nanotube through commonly named n- 
stacking interaction that we interpret here as weak and van der Waals interac- 
tions. Consequently, they are confined in the nanotube leading to different elec- 
tronic properties. As a preliminary step, we firstly present some results on the 
isolated coronene-coronene dimer interaction. In a second step, we will present 
the interaction of a stacked coronene chain encapsulated inside the (19,0) carbon 
nanotube. The chosen size of the nanotube corresponds to that experimentally ob- 
served by Okazaki and coworkers UJ, adequate to accomodate in a suitable way 
the coronene column, ii) In a second part, we will present the electronic properties 
determined by "many-body"-corrected DFT and the GW + Bethe-Salpeter calcu- 
lations. This framework will be used to determine the photoexcitation spectrum 
of an infinite clean coronene- stacking column, as well as the coronene column in- 
side the nanotube. As a complement, previously we also present a detailed study 
of the band-structure of this new hybrid system in order to interpret the obtained 
photoexcitation results. 

3.1. Coronene dimer 

We have first considered the interaction between two coronene molecules 
forming an isolated dimer (without the nanotube), in order to determine the "nat- 
ural" configuration of the dimer. All the calculations here have been performed 
using the LCAO-S* 2 + vdW formalism explained above. In Figure [I] (top panel) 
we have depicted the cohesive energy of the coronene dimer for various tilt an- 
gles as a function of the stacking distance between the molecules. The tilt angle 
is defined as a rotation angle of both molecules in parallel with respect to the 
y-axis, as it can be seen in the top sketch of Figure [TJ This specific angle will 
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Figure 2: (Color online). (Bottom panel) Cohesive energy per coronene molecule (inserted be- 
tween two neighboring coronenes) as a function of the intermolecular stacking distance of an 
infinite chain of stacked C24H 12 coronenes encapsulated inside the (19, 0) SWCNT. This cohesive 
energy is depicted for several tilt angles, considering a rotation angle of 30 degrees between the 
coronenes. The most stable structure (sketched in top panel) is obtained for a 35 degrees tilt an- 
gle, at an intermolecular stacking distance of around 4 A. In bottom panel, total energy minimum 
region is shown as an enlarged-scale inset for a better visualization. 
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be considered with caution to define the coronene stacking- structure in the nan- 
otube. We can observe a minimum energy of about -1.2 eV for the whole dimer, 
at a distance of 3.3 A, and for a tilt angle of 20 degrees. This result is rather 
different than for the coronene tilt angle in the monoclinic crystal, of around 46 
degrees for a 3.4 A intermolecular distance [|43l . since obviously we only consider 
an isolated coronene dimer. In a second calculation, we have determined, within 
the same formalism, the coronene dimer cohesive energy, considering the rota- 
tion of one molecule with respect to the other and with respect to the dimer-axis 
(z-axis), at tilt angle equal to 0. In this case, the minimum energy is found for 
a rotation angle of 30 degrees, and for a distance of 3.1 A, as represented in the 
middle panel of Figure [T] The geometrical configuration of the most stable situa- 
tion is also represented in the middle sketch of Figure [T] This result is confirmed 
by previous studies, where the specific 30 degrees configuration is denominated 
as "twisted-sandwich" structure [|45l . In a previous DFT study [46 ], by using a 
general-gradient B3LYP corrected functional, with no van der Waals correction, 
this structure is not found, however, to be the most stable one, showing the limi- 
tation of the standard DFT approach to describe these systems. In a recent work 
by our group, this essential role of the van der Waals interaction for the accurate 
determination of the geometrical structure has already been demonstrated for the 
Ceo fullerene confined in a carbon nanotube [44J. From a physical point of view, 
this stable situation corresponds to an A£>-stacking configuration, similarly to the 
graphene case, whilst in the previous calculation, at zero tilt angle, it corresponds 
to an AA-stacking, which is more repulsive. This behaviour is corroborated by 
the obtained cohesive energy, which in this case is of around -1.3 eV. 

Finally, in order to prepare the calculation of the stacked coronenes inside the 
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carbon nanotube, we have calculated the isolated coronene dimer cohesive energy 
for several tilt angles, comparing in each case the situation with rotation angles 
of and 30 degrees. Results are shown in bottom panel of Figure [TJ where we 
can clearly appreciate that the most stable configuration corresponds to a zero tilt 
angle with 30 degrees of rotation angle, obtaining the same stacking distance and 
energy as in the previous calculation. This seems to point out that the stacking 
of the coronenes inside the carbon nanotube will correspond to a deformation of 
this stable situation, due to the confinement within the nanotube walls. There- 
fore, in the next subsection, we will focus on this particular "twisted-sandwich" 
configuration to determine the coronenes packing inside the carbon nanotube. 

3.2. Encapsulated coronene-stacking inside the carbon nanotube 

In this subsection we present calculations of the encapsulation of stacked 
coronene molecules inside a (19,0) SWCNT. For this purpose, we have consid- 
ered 3 coronenes inside the nanotube, as represented in the top sketch of Figure [2| 
and we focus on the cohesive energy of the coronene lying in the middle. Figure [2] 
shows the evolution of this cohesive energy as a function of the inter-coronene 
stacking distance for various tilt angles, and considering a rotation of 30 degrees 
with respect to the tube axis from one coronene to the other. This rotation angle 
has been chosen according to the previous results obtained for the single dimer 
calculation. The most stable structure is obtained for a tilt angle of 35 degrees, at 
an inter-coronene distance of 4 A. The corresponding cohesive energy is around 
of 9.7 eV per coronene, with an excellent agreement with previous calculations 
on Ceo encapsulated inside a carbon nanotube [fT8l . This result is also in a good 
agreement with the experimental result giving a tilt angle of ~23 degrees, as well 
as with the corresponding DFT calculations giving a tilt angle of 35 degrees and 
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an inter-coronene distance of 4.2 A [1J. Nevertheless, as indicated above, the 
coronene stacking structure is rather different from the stability required by the 
system with a 30 degrees rotation of each coronene with respect to its neighbor- 
ing molecules. In order to check the validity of this result, we have also calcu- 
lated the same structure within our formalism, with no rotation angle between the 
coronenes. In this case we were not able to determine a stable configuration in 
the system, since the lowest energy is obtained for a 45 degrees tilt angle. We 
guess that in such a case, the most stable configuration should be obtained for a 
tilt angle approaching to 90 degrees, which is equivalent to have no stacking be- 
tween coronene molecules, but still forming chains. Of course, this result does not 
correspond at all to the Physics observed experimentally, meaning that the stable 
configuration requires a rotation of each coronene with respect to the other. There- 
fore, our study reveals a major aspect of coronene stacking encapsulated in carbon 
nanotubes, which is the 30 degrees rotation angle of each coronene with respect 
to the others. This aspect is even more important, since it has not been observed 
in this experiment. In addition to the calculation with a 30 degrees rotation, we 
represent in Figure[2]the cohesive energy of the same system with no rotation, and 
for a tilt angle of 45 degrees, which is the lowest energy we obtained in this case. 
It is noticeable the existing moderately high energy difference of around 1 eV per 
molecule with respect to our total minimum, taking into account the 30 degrees 
rotation angle, which validates our model. 

The discrepancy between our equilibrium structure and the one obtained in Q], 
especially regarding the 30 degrees rotation angle resulting in our calculations, 
may be explained in terms of two different factors: i) first, our treatment takes 
into account weak and van der Waals interaction, which was not the case in the 
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previous study, based on standard DFT calculations, and ii) second, in [1J, the 
authors have maintained a constant distance of 4.2 A between the coronenes in 
the nanotube, not allowing further lattice optimization. In the calculations we 
performed on this system, with no rotation angle, at a fixed distance of 4 A, we 
also determine that the 35 tilt angle yields the minimum energy of the system, 
eventhough this is not the absolute minimum (see Figure [2]). By varying the dis- 
tance, one can observe that in this configuration, increasing the tilt angle will 
increase the cohesive energy, probably until the coronenes stand all aligned inside 
the carbon nanotube. Consequently, our study goes beyond in the accurate char- 
acterization of this system, since we found that coronenes can only be stabilized 
by rotating each coronene with respect to its neighbors by 30 degrees. Since it 
has not been observed experimentally, due essentially to the low response of hy- 
drogen atoms in electron microscopy, we suggest future experimental work based 
on the observation by high-angle annular dark field scanning transmission elec- 
tron microscopy (HAADF-STEM) [|47ll48l of a similar system with some fluorine 
atoms substituting hydrogen atoms in the coronenes, as an example. One could 
consider stacked coronenes where some hydrogen atoms could be subtituted by 
Fluorine atoms, in order to be observed by microscopy techniques, in such a way 
to be able of discriminating (and categorically determining) the exact rotation an- 
gle between coronenes in their ground-state configuration. We have studied such 
system within DFT and observed that it is perfectly stable in the nanotube, with- 
out altering any of the structural properties. This procedure could be applied to a 
large amount of molecules stacked inside carbon nanotubes, since this rotation is 
probably a general behaviour in these systems. 
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4. "Many-body" Electronic Structure and Excitation Spectra 

In this section we will present the results of electronic structure calculations, 
as well as excitation spectrum, of an infinite coronenes chain encapsulated in the 
carbon nanotube, as to be compared with the infinite coronenes chain in its clean 
dilute-phase. For the isolated infinite coronenes chain (dilute phase) we have taken 
the most stable structure obtained for the dimer, periodically extended to mimic 
the experimental dimensionality, and with the optimum stacking distance obtained 
in the previous section. Nevertheless, the geometry taken for the simulation of 
the infinite coronenes chain encapsulated inside the nanotube is slightly different 
from the most stable one obtained from the exhaustive vdW analysis. Due to the 
high computational requirements of our GW+BSE formalism, for this system we 
have taken the ground-state structure for the coronenes encapsulated inside the 
nanotube, with no accounting the 30 degrees rotation angle. This rotation angle 
cannot be considered because the performance of our accurate formalism is just 
limited to one coronene per unit cell. The staking distance is that obtained in 
the geometrical analysis section. However, from our experience and regarding 
the next results, the rotation angle together with the electrostatic vdW interaction 
should not have a significant influence in the electronic structure scheme of the 
whole interacting system. 

4.1. Band- structure diagrams 

In Figure [3]<\ we show the Quasi-particle GW-corrected DFT band- structure 
diagram (referred to the Fermi energy level) for the isolated infinite (19, 0) zig-zag 
SWCNT for the most important symmetry line r — > X (along nanotube axis). The 
LDA+GW band diagram reveals a value for the band gap of e^ 9 ' ^ = 0.76 eV at 
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Figure 3: (Color online) Quasi-particle GW-corrected DFT band-structure diagrams (referred to 
the Fermi energy level) for: A) Isolated (19, 0) SWCNT. HOMO and LUMO levels of the self- 
assembled interacting C24H12 coronenes are also shown superimposed onto the nanotube band 
diagram (solid blue and red lines, respectively); and B) (19,0) SWCNT@C24H 12 interacting sys- 
tem, where HOMO and LUMO coronene states have been visually emphasized within the band 
diagram (solid blue and red lines, respectively), along the most relevant symmetry line from V — > 
X (along nanotube axis). Additionally, 3D orbital isodensities, all of them with a value of 2 x 10~ 5 
e~/A 3 , are shown for the most representative states of the systems. 
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the T fc-point, to be directly compared with those shown in previous literature of 
0.58 and 0.62 eV obtained by less accurate methods: i.e., third-order nonlinear 
optical response [49] and standard DFT by using a GGA-B3LYP hybrid potential 
(501, respectively. To our knowledge, no experimental band-gap value is available 
for this nanotube. As known, the obtained band gap manifests the narrow band- 
gap semiconducting character of the (19, 0) zig-zag SWCNT (for more details 
see |]49l [50l and references therein). In Figure [3j^., we also show superimposed 
on the band diagram the HOMO and LUMO levels of the self-assembled stacked 
C24H12 coronenes in their clean diluted phase (solid blue and red lines, respec- 
tively), keeping the same geometry and stacking periodicity they exhibit inside the 
nanotube fix. The obtained HOMO-LUMO gap has a value of e [ g 9 ' 0) = 2.92 eV, 
where the HOMO and LUMO levels (referred to the Fermi energy level) lie inside 
the unoccupied and occupied bands of the isolated tube, respectively, before the 
interaction between both systems is produced. For completion, 3D orbital HOMO 
and LUMO isodensities (both with a value of 2 x 10~ 5 e~/A 3 ) are depicted in Fig- 
ure^ The 3D isodensity profiles for the self-assembled stacked C 2 4Hi 2 coronenes 
visibly show a higher localization of the orbital HOMO charge around the double 
bonds inside the coronenes, whilst it is appreciable how the local symmetry of the 
LUMO completely changes with respect to the latest, making an optical transition 
between them permitted through the dispersion between the hybridized 7r-orbitals 
of neighboring coronenes, as explained in detail the spectroscopic analysis of ll35l . 
This analysis shall be extended in the excitation spectrum section. 

Figure [3jB shows the Quasi-particle GW-corrected DFT band- structure dia- 
gram (referred to the Fermi energy level) for the (19, 0) SWCNT@C 24 H 12 inter- 
acting system, where HOMO and LUMO coronene states have been visually em- 
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phasized within the band diagram (solid blue and red lines, respectively), along 
the most relevant symmetry line from r — > X (along nanotube axis). As an inter- 
esting result of the interaction between both subsystems, the LDA+GW band dia- 
gram reveals a global dilatation of all the bands around the Fermi energy, where a 
compression of occupied and empty bands is produced opening the original semi- 
conducting band gap of the tube up to a value of e^ 9 ' ^ = 3.4 eV at T point. Now, 
the interacting HOMO and LUMO levels of the self-assembled stacked C 2 4H 12 
coronenes (again solid blue and red lines) appear in the band diagram, but still, as 
expected, with an almost inexistent dispersion with the k points, mainly due to the 
low vdW interaction between them, and also with the wall of the tube. Neverthe- 
less, this apparently weak interaction of the internal coronenes with the tube is cer- 
tainly responsible of the nanotube gap-opening and the compression of the states 
in the global band- structure diagram. It is noticeable how the coronene-phase 
LUMO state remains located just like before the interaction is produced with re- 
spect to the Fermi energy, but, after the interaction, it interestingly overlaps in 
energy with the new nanotube conduction band. However, no apparent hybridiza- 
tion is produced between both states, each of them keeping its original coronene 
/ nanotube character (justified below). Meanwhile, the coronene-phase HOMO 
has moved up towards the Fermi energy level, by globally reducing their origi- 
nal HOMO-LUMO gap up to a value of e c ° r = 2.52 eV. Other interesting effect 
we can emphasize is that apparently the ^-dispersion profile of the nanotube bands 
seems not to show any variation before and after the interaction is produced, which 
manifest the robustness of the nanotube electronic structure, and the weak chem- 
ical interaction between the coronenes and the nanotube. Thus, summarizing, the 
interaction between the stacked coronenes inside the (19,0) nanotube produces 
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a three-folded effect: i) an opening of the nanotube band gap, losing the semi- 
conducting character; ii) a compression of the occupied and unoccupied bands 
around the Fermi level; and iii) a closing of the coronene-phase HOMO-LUMO 
gap, showing a noticeable hybridization between the coronene-phase LUMO state 
and the former nanotube conduction band. 

Similarly to the previous case, 3D isodensity profiles for the SWCNT@C 24 Hi 2 
system are also shown for the most relevant bands / states: former nanotube va- 
lence and conduction bands, and HOMO and LUMO states of the coronene-phase 
(all of them with a value of 2 x 10~ 5 e~/A 3 ). As appreciable in Figure [3} the 
3D isodensities associated to the nanotube valence and conduction band after 
the interaction is produced reveal a location of the orbital charge completely sur- 
rounding the tube. Regarding the same value taken in all cases for the depiction 
of the 3D isosurfaces, this justifies our previous statement, when we concluded 
no hybridization at all between the valence or conduction bands with any of the 
coronene-phase states; in particular, between the new nanotube conduction band 
and the coronene LUMO state which are located (not hybridized) at the same 
energy with respect to the Fermi level. Once again, 3D coronene-phase HOMO 
and LUMO orbital isodensities show a higher localization of the orbital HOMO 
charge around the double bonds inside the coronenes, with no variations with re- 
spect to their isolated diluted case, which indicated that the interaction between 
the coronenes and the nanotube is weak. In the same way, orbital LUMO iso- 
density profile keeps unaltered as well its symmetry after the interaction between 
the nanotube and the coronenes is produced. Giving the similar symmetry of the 
HOMO and LUMO states than for the non-interacting case, this interpretation for 
the interacting (19, 0) SWCNT@C 24 H 12 system still makes an optical transition 
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between them permitted. Nevertheless, the optical transition will be produced, 
in this case, at slightly lower excitation energy, since the HOMO-LUMO gap has 
been reduced by around ~ 0.4 eV. This analysis shall be extended in the excitation 
spectrum section. 

4.2. Excitation spectra 

Figure [4] shows the GW-BSE photoabsorption cross sections (calculated as 
explained above) as a function of the excitation energy for both the diluted self- 
assembled C 2 4Hi2 coronenes phase (solid blue line), and the (19, 0) SWCNT@C 24 Hi 2 
interacting system (solid green line) up to 6 eV. The finite width of the excitation 
peaks in an experimental spectrum - linked to the accessible resolution - is mostly 
determined by the temperature. In our calculations, on the other hand, each exci- 
tation peak in the spectrum has been broadened by a Lorentzian profile, where the 
broadening parameter r] is set equal to 0.2 eV, a value commonly used in order to 
mimic the experimental resolution [|42l . 

Regarding the figure, no significant photoabsorption signal is observed up to 
around 2 eV for both systems; nevertheless, the low-energy absorption regime (up 
to 2 eV) is slightly more pronounced for the SWCNT@C 24 H 12 interacting sys- 
tem, which may be due to the (weak) interaction between the stacked coronenes 
and the nanotube, and the dispersion of extended states coming from the periodic 
tube. The spectra of the self-assembled C 24 H 12 coronenes phase and the (19, 0) 
SWCNT@C 24 H 12 interacting system up to 4 eV are characterized by a main peak 
centered at 2.97 and 2.49 eV, respectively. According to the band diagrams shown 
in Figure |3j this peak seems to correspond to a permitted optical transition (see 
the orbital-based analysis below) from the HOMO state to the LUMO state of the 
staked coronenes chain (isolated or confined inside the nanotube). The intensity 
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Figure 4: (Color online) GW-BSE photoabsorption cross sections (in arbitrary units) as a function 
of the excitation energy (in eV) for: i) diluted self-assembled C24H12 coronenes phase (solid 
blue line); and ii) (19, 0) SWCNT@C2.4H12 interacting system (solid green line). In order to 
mimic the experimental resolution [49] the excitation peaks have been broadened by Lorentzian 
functions. The corresponding geometrical structure is shown in Figure [2] Photograph of the 
fluorescence from coronenesQSWCNTs in SDBS micelle solution (inset: left) and coronene in 
n-hexane solution (inset: right) under UV lamp irradiation (lex=356 nm) (2. 
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of this pronounced peak is slightly higher for the case of the coronenes inside the 
nanotube, which means that the electronic transition probability from HOMO to 
LUMO states increases for the confined coronenes. Nevertheless, both spectra 
profiles are quite similar, except for a relative red- shift of about 0.48 eV of the 
confined coronenes with respect to their dilute phase. Besides, both also show a 
very similar photoab sorption profile for the high absorption energy range (from 
4.5 eV), arising from less bounded states weaker in energy. 

The experimentally observed red-shift between the excitation spectra, roughly 
introduced in the previous section, is responsible of the change of the optical light- 
color emission from blue to green (see inset in Figure [4]). Photons of selected 
energies responsible to excite electrons from the HOMO to the LUMO states in 
the stacked coronenes, decay subsequently by the emission of blue or green light. 
As explained, the interaction between the stacked coronenes inside the (19, 0) 
nanotube produces a compression of the occupied and unoccupied bands around 
the Fermi level, and a closing of the coronene -phase HOMO-LUMO gap from 
2.97 eV (corresponding to a blue-light emission) to 2.49 eV (corresponding to a 
green-light emission). The origin of this optical red-shift from one system to the 
other can be justified in base to a detailed orbital-based analysis, which allows 
us to properly identify and label the HOMO and LUMO coronene-phase states 
(see isosurfaces in Figure [3]), respectively, as e2 g and 62 U in D 6 h symmetry (in 
spectroscopic nomenclature), and thus the transition between these states is per- 
mitted, but only allowed through z-axis polarized light (with a 2u symmetry) - in 
this case along the longitudinal stacking axis -, while it is forbidden through x and 
v-axis polarized light (with e lu symmetry). This interesting fact is also confirmed 
directly by the GW-BSE photoabsorption results, in which the z-axis polariza- 
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tion gives the only contribution to the lowest energy peaks in spectra of Figure |4j 
Moreover, it can be easily checked that there are no other allowed transition be- 
tween the HOMO of the coronene-phase and the frontier orbitals, apart from the 
mentioned permitted e 2g — > e 2u transition (see more details about this analysis 
in ll35l and references therein). This transition is produced in both systems at dif- 
ferent HOMO-LUMO gap energies of the coronene-phase, depending on if they 
are in the clean dilute phase, or confined inside the nanotube, which is manifested 
in the blue and green colored-light emissions observed in the experiment. 
Obviously, all these results as well as the van der Waals structures determined 
previously are sensitive to the particular choice of nanotube that we study here. In 
a future work we can think about a systematic study of coronene encapsulation in 
nanotube with respect to the chirality and the size of the host nanotube. 

5. Conclusions 

Following the experiment of Okazaki et al. Q, we have studied theoreti- 
cally the encapsulation of coronenes in a carbon nanotube, focusing on the (19, 0) 
SWCNT@C 2 4Hi2 system. In a first part, using a DFT-based formalims which 
takes into account weak and van der Waals interactions, we have studied the 
coronene dimer as a previous step for an accurate description of stacked coronenes 
encapsulated in the carbon nanotube. The major result that we obtain is a rotation 
of 30 degrees of each coronene with respect to its neighbors along the nanotube 
axis, leading to a tilt angle in the nanotube of 35 degrees and an intermolecular 
distance of 4 A. Actually this rotation between molecules, which has not been 
seen experimentally, is abolutely necessary to obtain the stability of the system. 
Whereas in previous DFT calculations (without any van der Waals interaction), 
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at fixed intermolecular stacking distance, the general configuration is very similar 
to our result, the 30 degrees rotation to ensure stability has not been shown be- 
fore. In a second step, starting from this determined configuration, we have used 
the GW formalism to accurately calculate the bandstructure of the whole (19, 0) 
SWCNT@C 24 Hi2 system and compared it with the isolated (19, 0) SWCNT. As a 
main result, we can clearly observe a compression of the coronene levels as well 
as an opening of the SWCNT gap, due to the mutual interaction between the two 
systems. Moreover, a precise determination of the excitation spectrum shows a 
clear red- shift from 2.97 eV in the coronene gas phase to 2.49 eV for the encap- 
sulated system. This result is in excellent agreement with the experimental results 
by Okazaki et al. 0], confirming the transition of the photoemission spectrum 
from blue to green emitted light. It can be deduced also that even if the stacking 
of coronene is a requisite for the stability of the structure (i.e. through the 30 
degrees rotation angle and the consequent tilt angle), the drastic change in elec- 
tronic structure is purely due to the interaction with the nanotubes and the band 
compression. Finally, we suggest here an experiment to observe this fundamental 
rotation between coronenes in the nanotube (which is probably a general trend in 
such molecular systems), by substituing in each coronene some of the hydrogen 
atoms by Fluorine atoms for example, these last ones being more observable in 
transmission electron microscopy. 
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